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The post-Newtonian approximation is a method for solving Einstein's 
field equations for physical systems in which motions are slow com- 
pared to the speed of light and where gravitational fields are weak. 
Yet it has proven to be remarkably effective in describing certain 
strong-field, fast-motion systems, including binary pulsars contain- 
ing dense neutron stars and binary black hole systems inspiraling 
toward a final merger. The reasons for this effectiveness are largely 
unknown. When carried to high orders in the post-Newtonian se- 
quence, predictions for the gravitational-wave signal from inspiraling 
compact binaries will play a key role in gravitational-wave detection 
by laser-interferometric observatories. 

general relativity | black holes | gravitational radiation 

Introduction 

The first detection of gravitational radiation, anticipated to 
occur during tliis decade, will undoubtedly be a triumph of 
experimental physics. It will initiate a new kind of astronomy. 
But it will also shine a spotlight on the theory of general rel- 
ativity itself. Gravitational radiation is a crucial prediction 
of Einstein's theory; indeed it is a natural prediction, given 
that the theory is built on a foundation of Lorentz invariance, 
which carries with it the concept of a limiting speed for inter- 
actions. 

On the other hand, general relativity is a notoriously com- 
plicated, non-linear, tensorial theory of the gravitational field. 
Almost no physically useful exact solutions of the theory are 
known, and those that are known are endowed with such high 
degrees of symmetry that their realm of validity is limited. To 
be sure, two of these solutions, due to Schwarzschild and Kerr, 
have proven to be of enormous importance, describing as they 
do the spacetime of isolated black holes, now widely accepted 
as being ubiquitous throughout the universe. Part of the 
Schwarzschild solution also describes the exterior geometry 
of any static spherical star or planet. Another extremely use- 
ful, albeit special solution is the Friedman-Robertson- Walker 
metric of the standard model of homogeneous and isotropic 
big-bang cosmology. 

But by its very nature, gravitational radiation involves 
spacetimes that are highly non-symmetrical and highly dy- 
namical. No exact solution of Einstein's equations is known 
that describes the emission and propagation of gravitational 
waves from a source, and the reaction of the source to the 
emission of those waves. 

As a result, most of our understanding of gravitational ra- 
diation has come from approximations to Einstein's equations. 
One class of approximations assumes that the gravitational 
fields in and around the source are suitably weak (the fields 
of the propagating waves weaken progressively as they leave 
the source), and that the motions within the source are suit- 
ably slow compared to the speed of light. This class includes 
schemes known as post-Newtonian theory, which will be the 
main subject of this paper, and a related scheme known as 
post-Minkowki theory. The underlying idea is to treat space- 



time as being that of flat Minkowski spacetime as the zeroth 
approximation, and to modify it by successive corrections. 

Another class of approximations takes a known exact so- 
lution of Einstein's equations, such as the black hole solutions 
of Schwarzschild or Kerr, and introduces small perturbations 
of those spacetimes, induced, for example, by a small par- 
ticle orbiting the hole. In cosmology, perturbations of the 
Friedman-Robertson- Walker solution permit treatment of the 
growth of large-scale structure in the universe and fluctuations 
in the cosmic background radiation. 

A rather different approach to solving Einstein's equations 
is "numerical relativity" , which endeavors to formulate and 
solve the exact equations to a precision limited only by avail- 
able computer resources for highly dynamical, highly asym- 
metrical situations using accurate and robust numerical com- 
putations. In recent years, numerical relativity has signifi- 
cantly enhanced our understanding of colliding black holes 
and neutron stars and the associated emission of gravitational 
radiation. 

Just as exact solutions of Einstein's equations have limited 
realms of validity, so too do approximation schemes. Specifi- 
cally, the post-Newtonian (FN) approximation is formally lim- 
ited to weak gravitational fields and slow motions. Yet recent 
experience has shown that the post-Newtonian approxima- 
tion is "unreasonably effective" in describing physical systems 
whose characteristics either lie outside the technical realm of 
validity of the approximation or push up against the boundary 
of that realm. 

The use of the term "unreasonably effective" , and indeed 
the title of this paper, have been shamelessly appropriated 
from a famous 1960 paper by Eugene Wigner, entitled "On 
the unreasonable effectiveness of mathematics in the physical 
sciences" 1:. In that paper, Wigner states "... the enormous 
usefulness of mathematics in the natural sciences is something 
bordering on the mysterious and . . . there is no rational expla- 
nation for it". Einstein presented similar ideas in an address 
to the Frussian Academy of Sciences in 1921. 

While the considerations of Einstein and Wigner deal with 
deep questions of the nature of knowledge of the physical 
world and how it is acquired and assembled, the topic of this 
paper is much more narrow. Our purpose will be to illus- 
trate the various ways in which the post-Newtonian approx- 
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imation has proven its extraordinary effectiveness in gravita- 
tional physics. Nevertheless it is no less mysterious: we have 
no good understanding of why this approximation to general 
relativity should be so effective. 

Because most of these applications of post-Newtonian the- 
ory involve the motions of gravitating bodies in their mutual 
gravitational fields and the emission of gravitational radiation, 
we will begin with a history of the "problem of motion" , one 
of the central challenges in the development of general rela- 
tivity. The problem has at times been contentious. We will 
then describe briefly the nature and structure of the post- 
Newtonian approximation, and will review some of the recent 
developments in the subject, particularly those involving go- 
ing to very high orders in the approximation. We will describe 
the usefulness of the PN approximation for characterizing al- 
ternative theories of gravity, and as a tool in experimental 
gravitation. Then we show how the PN approximation effec- 
tively and accurately accounts for the observed behavior of 
binary pulsar systems, despite the presence of strong gravi- 
tational fields in the interiors of the orbiting neutron stars. 
Finally, we will describe how the PN approximation, when 
carried to high enough orders, effectively describes the mo- 
tion of and gravitational radiation from inspiralling binary 
black hole systems, well into the strong-field region, and how 
its predictions merge smoothly with those from numerical rel- 
ativity. We will conclude by returning to the mystery of this 
effectiveness. 



Motion and radiation in general relativity: A history 

At the most primitive level, the problem of motion in general 
relativity is relatively straightforward, and was an integral 
part of the theory as proposed by Einsteirj^ A "test parti- 
cle", that is a particle whose mass is sufficiently small that 
its own contribution to the curvature of spacetime can be ig- 
nored, moves on a geodesic of the curved spacetime in which it 
finds itself. Underlying this concept is the "weak equivalence 
principle" , which states that the acceleration of a suitably 
small body in an external gravitational field is independent 
of its internal structure or composition, a principle verified 
by modern experiments to parts in 10^^. Using the geodesic 
equation and an approximate solution for the spacetime met- 
ric around the Sun, Einstein was able in 1915 to obtain the 
relativistic contribution to the perihelion advance of Mercury, 
in agreement with observations. 

The first attempts to treat the motion of multiple bodies, 
each with a finite mass, were made in the period 1916-1917 by 
de Sitter, Lorentz and Droste [S] S] . They derived the metric 
and equations of motion for a system of A'^ bodies, in what 
today would be called the first post-Newtonian approxima- 
tion of general relativity (de Sitter's equations turned out to 
contain some important errors). In 1916, Einstein took the 
first crack at a study of gravitational radiation, deriving the 
energy emitted by a body such as a rotating rod or dumbbell, 
held together by non-gravitational forces [Q. He made some 
unjustified assumptions as well as a trivial numerical error 
(later corrected by Eddington [6J), but the underlying con- 
clusion that dynamical systems would radiate gravitational 
waves was correct. 

The next significant advance in the problem of motion 
came 20 years later. In 1938, Einstein, Infeld and Hoff- 
man published the now legendary "EIH" paper, a calculation 
of the A-body equations of motion using only the vacuum 
field equations of general relativity [?]• They treated each 
body in the system as a spherically symmetric object whose 
nearby vacuum exterior geometry approximated that of the 
Schwarzschild metric of a static spherical star. They then 



solved the vacuum field equations for the metric between each 
body in the system in a weak field, slow-motion approxima- 
tion. Then, using a primitive version of what today would be 
called "matched asymptotic expansions" they showed that, in 
order for the nearby metric of each body to match smoothly 
to the interbody metric at each order in the expansion, cer- 
tain conditions on the motion of each body had to be met. 
Together, these conditions turned out to be equivalent to the 
Droste-Lorentz A-body equations of motion. Remarkably, the 
internal structure of each body was irrelevant, apart from the 
requirement that its nearby field be approximately spherically 
symmetric. 

Around the same time, there occurred an unusual detour 
in the problem of motion. Using equations of motion based 
on de Sitter's paper, specialized to two bodies, Levi-Civita [8] 
showed that the center of mass of a binary star system would 
suffer an acceleration in the direction of the pericenter of the 
orbit, in an amount proportional to the difference between 
the two masses, and to the eccentricity of the orbit. Such an 
effect would be a violation of the conservation of momentum 
for isolated systems caused by relativistic gravitational effects. 
Levi-Civita even went so far as to suggest looking for this ef- 
fect in selected nearby close binary star systems. However, 
Eddington and Clark 9 quickly pointed out that Levi-Civita 
had based his calculations on de Sitter's flawed work; when 
correct two-body equations of motion were used, the effect 
vanished, and momentum conservation was upheld. Robert- 
son confirmed this using the EIH equations of motion |10j . 
Ironically, the acceleration of the center of mass of a binary 
system is today a subject of great theoretical and astrophysi- 
cal interest, albeit for a very different reason. We will return 
to this subject later. 

Roughly 20 more years would pass before another major 
attack on the problem of motion. This period was the con- 
tinuation of a time of relative dormancy for the entire sub- 
ject of general relativity that lasted from the 1920s until the 
1960s. This dormancy resulted in part from the lack of exper- 
imental or observational relevance for the theory, in part from 
the perceived complexity of the theory, and in part from the 
emergence of new fields of physics such as nuclear and particle 
physics in the middle part of the 20th century. 

But in the middle 1960s, when a revival of general rela- 
tivity was in its early phase, Fock in the USSR and Chan- 
drasekhar in the US independently developed and system- 
atized the post-Newtonian approximation in a form that laid 
the foundation for modern post-Newtonian theory |lll 112] . 
They developed a full post-Newtonian hydrodynamics, with 
the ability to treat realistic, self-gravitating bodies of fluid, 
such as stars and planets. In the suitable limit of "point" par- 
ticles, or bodies whose size is small enough compared to the 
interbody separations that flnite-size effects such as spin and 
tidal interactions can be ignored, their equations of motion 
could be shown to be equivalent to the EIH and the Droste- 
Lorentz equations of motion. Chandrasekhar and his students 
also began extending the theory to higher orders in the post- 
Newtonian approximation. 

An important byproduct of the Fock-Chandrasekhar work 
was the discovery by Nordtvedt 13 that, in theories of gravity 
alternative to general relativity, the motion of self-gravitating 
bodies could depend on their internal structure, in contrast 
to general relativity, where a body's internal structure is ir- 
relevant to its motion (barring tidal or spin effects). As a 
result of this "Nordtvedt effect", the Earth and the Moon 
could fall toward the Sun with slightly different accelerations. 
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because of the small difference in their internal gravitational 
binding energy per unit mass. Nordtvedt's discovery led to 
an important new test of general relativity using Lunar laser 
ranging, and to the development by Nordtvedt and by Will of 
the parametrized post-Newtonian (PPN) framework for treat- 
ing alternative theories and experimental tests. This will be 
the subject of a later section. 

The next important period in the history of the problem of 
motion was 1974 -1979, initiated by the 1974 discovery of the 
binary pulsar PSR 1913-1-16 by Hulse and Taylor [2]. Around 
the same time there occurred the first serious attempt to cal- 
culate the head-on collision of two black holes using purely 
numerical solutions of Einstein's equations, by Smarr and col- 
laborators ^15, , building on the pioneering work by Hahn and 
Lindquist [16] 

The binary pulsar consists of two neutron stars, one an 
active pulsar detectable by radio telescopes, the other very 
likely an old, inactive pulsar. Each neutron star has a mass 
of around 1.4 solar masses. The orbit of the system was seen 
immediately to be quite relativistic, with an orbital period of 
only eight hours, and a mean orbital speed of 200 km/s, some 
four times faster than Mercury in its orbit. Within weeks of its 
discovery, numerous authors pointed out that PSR 1913-1-16 
would be an important new testing ground for general rela- 
tivity. In particular, it could provide for the first time a test 
of the effects of the emission of gravitational radiation on the 
orbit of the system. 

However, the discovery revealed an ugly truth about the 
"problem of motion". As Ehlers et al. pointed out in an in- 
fiuential 1976 paper J17', the general relativistic problem of 
motion and radiation was full of holes large enough to drive 
trucks through. They pointed out that most treatments of 
the problem used "delta functions" as a way to approximate 
the bodies in the system as point masses. As a consequence, 
the "self-field" , the gravitational field of the body evaluated 
at its own location, becomes infinite. While this is not a ma- 
jor issue in Newtonian gravity or classical electrodynamics, 
the non-linear nature of general relativity requires that this 
infinite self-field contribute to gravity. In the past, such infini- 
ties had been simply swept under the rug. Similarly, because 
gravitational energy itself produces gravity it thus acts as a 
source throughout spacetime. This means that, when calcu- 
lating radiative fields, integrals for the multipole moments of 
the source that are so useful in treating radiation begin to di- 
verge. These divergent integrals had also been routinely swept 
under the rug. Ehlers et al. further pointed out that the true 
boundary condition for any problem involving radiation by 
an isolated system should be one of "no incoming radiation" 
from the past. Connecting this boundary condition with the 
routine use of retarded solutions of wave equations was not a 
trivial matter in general relativity. Finally they pointed out 
that there was no evidence that the post-Newtonian approx- 
imation, so central to the problem of motion, was a conver- 
gent or even asymptotic sequence. Nor had the approximation 
been carried out to high enough order to make credible error 
estimates. 

During this time, some authors even argued that the 
"quadrupole formula" for the gravitational energy emitted by 
a system (see below), while correct for a rotating dumbell as 
calculated by Einstein, was actually wrong for a binary sys- 
tem moving under its own gravity. The discovery in 1979 
that the rate of decay of the orbit of the binary pulsar was in 
agreement with the standard quadrupole formula made some 
of these arguments moot. Yet the question raised by Ehlers 



et al. was still relevant: is the quadrupole formula for binary 
systems an actual prediction of general relativity? 

Motivated by the Ehlers et al. critique, numerous workers 
began to address the holes in the problem of motion, and by 
the late 1990s most of the criticisms had been answered, par- 
ticularly those related to divergences. The one question that 
remains open is the nature of the post-Newtonian sequence; 
we still do not know if it converges, diverges or is asymptotic. 
Despite this, it has proven to be remarkably effective. 

The final important development in this history of the 
problem of motion was the proposal to build large-scale laser 
interferometric gravitational-wave observatories, both on the 
ground and in space [181 [19 , 20 . It was quickly realized that, 
in order to maximize the likelihood of detecting the leading 
candidate source of gravitational waves - the final inspiral and 
merger of binary systems of neutron stars or black holes, ex- 
tremely accurate theoretical predictions of the gravitational 
waveform signal would be needed. This meant that calcula- 
tions of the equations of motion and gravitational waves from 
binary systems would have to be carried out to many orders 
beyond the simple post-Newtonian approximation. The com- 
pletion of this ambitious program by many groups worldwide 
has led to the discovery of the unreasonable effectiveness of 
post-Newtonian theory in the extreme realm of the merger of 
compact astronomical bodies. 

The post-Newtonian approximation 

The post-Newtonian approximation is based on the assump- 
tion that gravitational fields inside and around bodies are 
weak and that characteristic motions of matter are slow com- 
pared to the speed of light. This means that one can charac- 
terize the system in question by a small parameter e, where 

{v/cf ^ GM/rc^ ^ p/pc , [1] 

where v, M and r denote the characteristic velocity, mass, 
and size or separation within the system; p and p are the 
characteristic pressure and density within the bodies; G and 
c are Newton's gravitational constant and the speed of light, 
respectively. 

One then incorporates this approximation into methods 
for solving Einstein's equations. Those equations, G^^ = 
87r(G/c*)T'^y are elegant and deceptively simple, showing how 
geometry (in the form of the Einstein tensor G^ji^, which is a 
function of spacetime curvature) is generated by matter (in 
the form of the material energy-momentum tensor T^^). How- 
ever, this is not the most useful form for actual calculations. 
For post-Newtonian calculations, a far more useful form is the 
set of so-called "relaxed" Einstein equations: 

Uh°'^ = -167r(G/c*)T"'', [2] 

where □ = —d^/dictY' -\- is the fiat-spacetime wave oper- 
ator, h°'^ is a "gravitational tensor potential" related to the 
deviation of the spacetime metric g^fi from its flat-spacetime 
Minkowski form 7?„/3 by the formula h"'^ = - {-gf'^g"'^, 
where g is the determinant of gap, and where a particular 
coordinate system has been specified by the de Donder or 
harmonic gauge condition dh"^ /dx^ = (summation on re- 
peated indices is assumed, with x^ — ct). Because we as- 
sume that gravity is weak everywhere, the field h°'^ is "small" . 
This form of Einstein's equations bears a striking similarity 
to Maxwell's equations for the vector potential A" in Lorentz 
gauge: DA" = -47rJ", dA^/dx" = 0. The key difference is 
that the source on the right hand side of Eq. S is given by 
the "effective" energy-momentum pseudotensor, 

t"" = {-g)T°'^ + (167r)-'A"'^, [3] 
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where A"'' is the non-hnear "field" contribution given by 
terms quadratic (and higher) in /i"'' and its derivatives 
(see [21], Eqs. (20.20, 20.21) for formulae). In general rel- 
ativity, the gravitational field itself is a source of gravity, a 
refiection of the nonlinearity of Einstein's equations, and in 
contrast to the linearity of Maxwell's equations. 

Equation (|2| is an exact restatement of Einstein's equa- 
tions, and depends only on the assumption that spacetime can 
be covered by harmonic coordinates. It is called "relaxed" 
because it can be solved formally as a functional of source 
variables without specifying the motion of the source, in the 
form 



/i"''(t,x) = 



4G 



'(t - |x - x'|/c,x'; 

X — x' 



j3 ' 

a X , 



[4] 



where the integration is over the past fiat-spacetime null cone 
C of the field point (t,x). The motion of the source is then 
determined either by the equation dr"^ /dx^ = (which fol- 
lows from the harmonic gauge condition), or from the usual 
covariant equation of motion V pT"^ = 0, where V/3 denotes a 
covariant derivative. This formal solution can then be iterated 
in a weak-field (||/i°''^|| ^ 1) approximation. One begins by 
substituting h^^ = into the source t"^ in Eq. Q, and solv- 
ing for the first iterate /i" , and then repeating the procedure 
sufficiently many times to achieve a solution of the desired ac- 
curacy. This procedure is often called post-Minkowski theory. 
If one further imposes the slow motion {v <C c) assumption, 
one obtains post-Newtonian theory. For example, to obtain 
the equations of motion at first post-Newtonian (IPN) or- 
der, two iterations are needed (i.e. /ij^ must be calculated). 
To obtain the leading gravitational waveform and energy fiux 
far from a binary system, two iterations are needed, while 
to obtain the leading contributions to gravitational radiation 
damping of the source, three iterations are necessary (see [55] 
for a discussion). 

However, because the source t°'^ contains h"'^ itself, it is 
not confined to a compact region, but extends over all space- 
time. As a result, as emphasized by Ehlers et al. [T7], there 
is a danger that the integrals involved in various expansions 
of the solutions for h"^ will diverge or be ill-defined. Nu- 
merous approaches were developed to handle this difficulty. 
The "post-Minkowski" method of Blanchet, Damour and Iyer 
solves Einstein's equations by two different techniques, one 
in the near zone (within one gravitational wavelength of the 
source) and one in the far zone, and uses the method of sin- 
gular asymptotic matching to join the solutions in an overlap 
region. The method provides a natural "regularization" tech- 
nique to control potentially divergent integrals (see [23| for a 
thorough review of this method). The "Direct Integration of 
the Relaxed Einstein Equations" (DIRE)approach of Will, 
Wiseman and Pati (see [24|) retains Eq. Q as the global so- 
lution, but splits the integration into one over the near zone 
and another over the far zone, and uses different integration 
variables to carry out the explicit integrals over the two zones. 
In the DIRE method, all integrals are finite and convergent. 

The problem of "delta-functions" was handled in a variety 
of ways. One was to import from quantum field theory a set of 
powerful techniques called "dimensional regularization", that, 
subject to some benign assumptions about analyticity, could 
be used to control the potential infinities order by order in the 
PN expansion. Another adapted the EIH matching method to 
higher-order PN calculations. A third approach treated each 
body as a real, nearly spherical fiuid ball, and sorted contri- 
butions depending on how they scaled with the size of the 
ball. In every case where different methods made a predic- 



tion about the equations of motion or the gravitational-wave 
signal, they were in complete agreement. 

Among the results of these approaches are formulae for the 
equations of motion and gravitational waveform of binary sys- 
tems of compact objects, carried out to high orders in a PN 
expansion. Here we shall only state a few formulae for the 
purpose of illustration. For example, the relative two-body 
equation of motion has the form 

dv Gmi ^ 1 ^ 1 ^ 1 ^ 

— n H — jAipN H — -j;A2PN H — g-A2.5PN 



dt 



1 . 1 . 

+ -J-A3PN H — ^As.spN + . 



where m = mi -I- 7712 is the total mass, r = |xi — X2I, 
V = vi — V2, and n — (xi — X2)/r. The notation A„pn 
indicates that the term is 0(e") relative to the leading New- 
tonian term — n. Explicit and unambiguous formulae for non- 
spinning bodies through 3.5PN order have been calculated by 
various authors, and a number of spin-orbit and spin-spin con- 
tributions have been obtained (see [23] for a review). Here we 
quote only the first PN correction and the leading radiation- 
reaction term at 2.5PN order: 



AiPN = <^ (4 + 2r?) 



(1 + 3ri)v'^ + -rir" 



-H(4-2r?)rv, 



8 Gm 

A2.5PN = —rr^V 

15 r 



Qv" + 17 



Gm \ 



[7] 



where r) = mim2/(mi -I- 7712)^ and r = dr/dt. These 
terms are sufficient to analyse the orbit and evolution of bi- 
nary pulsars. For example, the IPN terms are responsible 
for the periastron advance of an eccentric orbit, given by 
dbj/dt — Q-K fhGm/[a(? {\ — e^)], where a and e are the semi- 
major axis and eccentricity, respectively, of the orbit, and /b 
is the orbital frequency, given to the needed order by Kepler's 
third law 27r/b = {Gm/a^Y''^. 

Another product is a formula for the gravitational field far 
from the system, whose spatial components K'-' (often called 
the gravitational "waveform") are sufficient to determine the 
signal detected by a laser interferometer, written schemati- 
cally in the form 



2Gm / 



c 



^O.SPN + 2 ^ 



;iPN ^ ^3 



,3 "^LSPN 



+ g4 Q2PN + ^5 Q 



2.5PN 1^ 



[81 



where R is the distance from the source, and the variables are 
to be evaluated at retarded time t — R/c. The leading term is 
the so-called quadrupole formula, given explicitly by 

h^'{t,^)^^r'{t^R/c), [9] 

where P-' is the quadrupole moment of the source, and over- 
dots denote time derivatives. For a binary system this leads 
to Q'-' = 2rj(v^v-' — Gmii'^n' /r). For binary systems, explicit 
formulae for the waveform through 3.5PN order have been 
derived (see |23| for a full review). 

Given the gravitational waveform, one can compute the 
rate at which energy is carried off by the radiation. The 
lowest-order quadrupole formula leads to the gravitational 
wave energy fiux 



E ■ 



8 2GW,^„ 2 



15 



llT^^ 



fioi 
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This has been extended to 3.5PN order beyond the quadrupole 
formula [53]. Formulae for fluxes of angular and linear mo- 
mentum can also be derived. The 2.5PN radiation-reaction 
terms in the equation of motion Q result in a decrease of 
the o rbit al energy at a rate that precisely balances the energy 
flux ( |10[ ) determined from the waveform. Averaged over one 
orbit, this results in a rate of increase of the binary's orbital 
frequency given by 



dt 



F{e) = (1 - 



-7/2 



, 73 2 37 4 
^+24^ +96^ 



[111 



3/5. 



where M is the so-called "chirp" mass, given by = r; 
Notice that by making precise measurements of the phase 
^{t) = 2n J* f{t')dt' of either the orbit or the gravitational 
waves (for which / — 2/b for the dominant component) as a 
function of the frequency, one in effect measures the "chirp" 
mass of the system. 



The parametrized post-Newtonian frameworl< and 
tests of general relativity 

The post-Newtonian approximation has been remarkably ef- 
fective as a tool for interpreting experimental tests of general 
relativity. This is because, in a broad class of alternative met- 
ric theories of gravity, it turns out that only the values of a 
set of numerical coefficients in the post-Newtonian expression 
for the spacetime metric vary from theory to theory. Thus 
one can encompass a wide range of alternative theories by 
simply introducing arbitrary parameters in place of the nu- 
merical coefflcients. This idea dates back to Eddington in 
1922, but the "parametrized post-Newtonian (PPN) frame- 
work" was fully developed by Nordtvedt and by Will in the 
period 1968-72 25, 26, 27 . The framework contains ten PPN 
parameters: 7, related to the amount of spatial curvature gen- 
erated by mass; /3, related to the degree of non-linearity in the 
gravitational fleld; ^, ai, 02, and 03, which determine whether 
the theory predicts that local gravitational experiments could 
yield results that depend on the location or velocity of the 
reference frame; and ^1, ^2, Cs ^^^^ Cii which describe whether 
the theory has appropriate momentum conservation laws. In 
general relativity, 7 = 1,13= 1, and the remaining parameters 
all vanish. For a complete exposition of the PPN framework 
see [28]. 

To illustrate the use of these PPN parameters in experi- 
mental tests, we cite the deflection of light by the Sun, an ex- 
periment that made Einstein an international celebrity when 
the sensational news of the Eddington-Crommelin eclipse 
measurements was relayed in November 1919 to a war-weary 
world. For a light ray which passes a distance d from the Sun, 
the deflection is given by 



A6' 



H-7\ 4GM 
2 I dc^ 



[12] 



where M is the mass of the Sun. The "1/2" part of the coeffi- 
cient can be derived by considering the Newtonian deflection 
of a particle passing by the Sun, in the limit where the parti- 
cle's velocity approaches c; this was first calculated indepen- 
dently by Henry Cavendish and Johann von Soldner around 
1800 (see, e.g., .29i). The second "7/2" part comes from the 
bending of "straight" lines near the Sun relative to lines far 
from the Sun, as a consequence of the curvature of space. A 



related effect called the Shapiro time delay, an excess delay in 
travel time for light signals passing by the Sun, also depends 
on the coefficient (1 + 7)/2. 

Measurements using visible light, made during solar 
eclipses, began with the 1919 measurements of Eddington and 
his colleagues, but never reached better than ten per cent 
precision, largely because of the logistical difficulties inher- 
ent in such measurements. High precision measurements were 
achieved using radio waves beginning in the late 1960s, culmi- 
nating in the use of Very Long Baseline Radio Interferometry 
(VLSI). A 2004 analysis of VLBI data on 541 quasars and 
compact radio galaxies distributed over the entire sky verified 
general relativity at the 0.02 percent level 30 . Shapiro time 
delay measurements began also in the late 1960s, by bounc- 
ing radar signals off Venus and Mercury or by tracking inter- 
planetary spacecraft; the most recent test used tracking data 
from the Cassini spacecraft while it was en route to Saturn, 
yielding a result at the 0.001 percent level [3T]. For a com- 
prehensive review of the current status of experimental tests 
of GR, see [31]. 

Other experimental bounds on the PPN parameters 
came from measurements of the perihelion-shift of Mercury, 
searches for the "Nordtvedt effect" in the Earth-Moon orbit 
using Lunar laser ranging, and a variety of geophysical and 
astronomical observations. All bounds were consistent with 
the predictions of general relativity, as summarized in Table 
1. 



Binary pulsars and the Strong Equivalence Principle 

Binary pulsars, such as the famous Hulse-Taylor system PSR 
1913+16 [13], illustrate the unreasonable effectiveness of the 
post-Newtonian approximation. 

Through precise timing of the pulsar "clock" , the orbital 
parameters of the system can be measured with exquisite pre- 
cision. These include non-relativistic "Keplerian" parameters, 
such as the orbital eccentricity e, and the orbital period Pb, as 
well as a set of relativistic, or "post-Keplerian" parameters. 
The latter parameters include the mean rate of advance of pe- 
riastron (duj/dt), the analogue of Mercury's perihelion shift; 
the effect of special relativistic time-dilation and the gravita- 
tional redshift on the observed phase or arrival time of pulses, 
resulting from the pulsar's orbital motion and the gravita- 
tional potential of its companion (represented by a parameter 
7'); the rate of decrease of the orbital period (dPt/dt), taken 
to be the result of gravitational radiation damping (apart from 
a small correction due to galactic differential rotation); and 
two parameters related to the Shapiro time delay of the pul- 
sar signal as it passes by the companion. According to general 
relativity, the post-Keplerian effects depend only on e and Pb, 
which are known, and on the two stellar masses, which are 
unknown a priori. By combining the observations of PSR 
1913+16 with the general relativity predictions for the first 
three post-Keplerian parameters, one obtains both a mea- 
surement of the two masses, and a test of the theory, since 
the system is overdetermined. The results are 

mi = 1.4414 ±O.OOO2Af0 ,7712 = 1.3867 ± 0.0002Mo , [13] 



= 1.0013 ± 0.0021 . 



[14] 



The concordance among the three constraints on the two 
masses is shown in Fig. [l] The accuracy in measuring the 
relativistic damping of the orbital period is now limited by 
uncertainties in our knowledge of the relative acceleration be- 
tween the solar system and the binary system as a result of 
galactic differential rotation. In the recently discovered "dou- 
ble pulsar" J 0737-3039, all five post-Keplerian parameters 
are measured, together with the mass ratio mi/m2 derived 
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directly from the ability to observe the motion of both pul- 
sars. All six constraints on the masses overlap one another on 
the mi — 7712 plane, again consistently with general relativity. 

However, there is something potentially wrong here. All 
the post-Keplerian effects discussed above are calculated using 
post-Newtonian theory. Yet the neutron stars that compose 
these systems have very strong internal gravity. This gravi- 
tational binding energy reduces the total mass of each body 
by 10 to 20 percent compared to the total rest mass of its 
constituent particles. By contrast, the orbital energy is only 
10"^ of the mass-energy of the system. Since general relativ- 
ity is a non-linear theory, surely there is some mixing between 
the strong internal gravity and the weak interbody gravity. So 
how can post-Newtonian theory possibly give valid predictions 
for such systems? 

The reason is a remarkable property of general relativity 
called the Strong Equivalence Principle (SEP). A consequence 
of this principle is that the internal structure of a body is "ef- 
faced" , so that the orbital motion and gravitational radiation 
emitted by a system of well separated bodies depend only on 
the total mass of each body, and not on its internal struc- 
ture, apart from standard tidal and spin-coupling effects. In 
other words, the motion of a normal star or a neutron star 
or a black hole depends on the body's total mass, and not on 
the strength of its internal gravitational fields. This behavior 
was already implicit in the work of Einstein, Infeld and Hoff- 
man, where only the exterior nearby field of each body was 
needed, and has been verified theoretically to at least second 
post-Newtonian order by more modern methods. 

By contrast, in alternative theories of gravity, SEP is not 
valid in general, and internal-structure effects can lead to sig- 
nificantly different behavior, such as the Nordtvedt effect, a 
possible difference in acceleration of the Earth compared to 
the Moon in the solar gravitational field, or such as the emis- 
sion of dipole gravitational radiation from systems of bodies 
with dissimilar internal structure. The close agreement of 
binary-pulsar data with the predictions of general relativity 
constitutes a kind of "null" test of the effacement of strong- 
field effects in that theory. It also constitutes a verification 
of the unreasonable effectiveness of post-Newtonian theory in 
this class of strong-gravity systems. 



Gravitational waves and inspiralling compact binaries 

Possibly the most remarkable example of the unreasonable ef- 
fectiveness of post-Newtonian theory is that of the inspiral 
and merger of binary systems of compact objects such as neu- 
tron stars and black holes. The decay of the orbit of a com- 
pact binary system through gravitational-wave emission will 
ultimately bring the two bodies together in a final and catas- 
trophic inspiral, followed by a merger and the likely formation 
of a terminal black hole. This process will emit a character- 
istic gravitational-wave signal with rising frequency and am- 
plitude (often called a "chirp") that should be detectable by 
the world-wide network of ground-based laser-interferometric 
observatories that is expected to be once again operational fol- 
lowing major upgrades, by 2015 |18lll9j . In the low-frequency 
end of the gravitational-wave spectrum, the merger of super- 
massive black holes in the centers of galaxies will be detectable 
from cosmological distances by the proposed space-based in- 
terferometer LISA, currently being planned for a launch after 
2020 [SOj, and by next-generation arrays of radio telescopes 
doing pulsar timing 33]. 

The most effective technique for detecting potential binary 
inspiral gravitational wave signals embedded in the noisy out- 
put of these interferometers is the method of matched filtering, 



whereby a theoretically generated gravitational-wave signal 
appropriate for a given source is cross correlated against the 
output of the detector. Since the noise is a random process, 
such a cross-correlation will yield a positive signature if there 
is a signal that precisely matches the template over the hun- 
dreds to thousands of cycles of signal that are expected to 
lie within the detectable band, even if the signal is formally 
weaker than the noise. With a bank of template waveforms 
that depend on the source parameters such as the two masses, 
spins, sky location, orbital eccentricity, etc, it will be possible 
both to detect signals and to measure the properties of the 
source [M]. To be most effective, this method requires very 
accurate theoretical templates. For the inspiral part of the 
signal, these templates have been calculated by many groups 
using post-Newtonian theory, with equations of motion (5| 
and gravitational waveforms Q calculated to 3.5PN order 
beyond the leading terms [23j . 

Eventually, however, the inspiral will reach the state where 
the orbital velocities are high and the gravitational fields are 
strong, so that the post- Newtonian approximation is no longer 
valid. Given our lack of knowledge of the convergence proper- 
ties of the approximation, it is not known a priori where this 
should occur. 

Meanwhile, it has become increasingly clear that the sig- 
nal from the final few inspiral orbits, from the merger of the 
two bodies, and even from the final vibrations of the newly 
formed black hole will make important contributions to the 
detectability of the waves by the interferometers. Luckily, 
many years after its primitive beginnings in the 1970s, numer- 
ical relativity finally reached a stage, following critical break- 
throughs in 2005 35 , where researchers could reliably and 
robustly simulate that final part of the inspiral process. 

It therefore came as a complete surprise, when gravita- 
tional waveforms from post-Newtonian theory were compared 
with numerical relativity waveforms for those final orbits, to 
discover that the agreement was unreasonably good. The am- 
plitudes and phases of the waves calculated by the two meth- 
ods agreed remarkably well cycle by cycle over many cycles, 
and this was in a regime where Gmjrc? ~ 0.2 and vjc ~ 0.4, 
where one had no right to expect the post-Newtonian approx- 
imation to be valid, even with many high-order correction 
terms in the PN formulae 36, 37, 38 . 

This unreasonable agreement was crucial, because it per- 
mitted the development of techniques for "stitching together" 
post-Newtonian and numerical relativity waveforms to obtain 
templates that are accurate and valid over the entire inspi- 
ral and merger process. Selecting the best stitching method 
involves taking into account the noise characteristics of the in- 
terferometers whose data is to be analysed, and to find stitch- 
ings that optimize all the data analysis protocols, such as 
false alarm thresholds, detection confidence criteria, and so 
on, that are part and parcel of all signal detection strategies. 
Because the post-Newtonian waveforms are analytic expan- 
sions, they can be resummed using such tricks as Pade ap- 
proximants to suggest alternative ways to match numerical 
results. This ongoing work involves a unique collaboration 
among post-Newtonian theorists, numerical relativists and in- 
terferometer data analysts (see, e.g. [39]), but it would have 
been moot had not post-Newtonian theory been so effective 
in overlapping with numerical relativity in the strong-gravity 
high-speed regime. 

Another example of the unreasonable effectiveness of post- 
Newtonian theory relates to the "kick" given to a black hole 
formed from the merger of two compact objects. In contrast 
to the erroneous 1937 claim by Levi-Civita, this is a real ef- 
fect. If a system emits gravitational waves anisotropically, 
then the waves carry linear momentum away in addition to 
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energy and angular momentum, and, by virtue of the over- 
all conservation of momentum, the source must recoil in the 
opposite direction. It turns out that the gravitational-wave 
recoil imparted to a final black hole could have important 
astrophysical consequences, especially for the mergers of su- 
permassive black holes, possibly ejecting the black hole from 
the host galaxy, or displacing it sufficiently from the center to 
cause interactions with surrounding gas or stars, thus gener- 
ating an electromagnetic counterpart signal '40'. 

Using formulae for the radiated momentum flux valid to 
2PN order, Blanchet et al. |4lj calculated the kick imparted 
to a black hole from the merger of two non-spinning black 
holes, as a function of their mass difi'erence (for equal masses, 
the effect vanishes by symmetry). The resulting kick velocity 
turned out to be in remarkable agreement with kicks deter- 
mined subsequently using numerical relativity, but only up to 
the point where the two black holes were about to merge. The 
agreement was surprising, since the dominant contribution to 
the recoil comes from a regime where the post-Newtonian ap- 
proximation should have failed. 

The numerical simulations also showed that, following the 
merger, there was a small "anti-kick", reducing the final re- 
coil velocity by around 30 percent Le Tiec et al. [^3] 
then used a hybrid calculation that combined formulae for the 
metric surrounding two closely spaced black holes accurate to 
2PN order as initial data, combined with the method of black 
hole perturbation theory, to study whether the anti-kick was 
produced by the linear momentum radiated during the ring- 
down phase of the vibrating final black hole. The analysis 
found that the combination of the 2PN inspiral kick plus the 
ringdown kick was in agreement with the kicks obtained by 
numerical relativity (Fig. [2|. 

A final instance of the effectiveness of post- Newtonian the- 
ory came from analyses of the initial configurations used in nu- 
merical relativity to study compact binary inspiral. Because 
gravitational radiation tends to circularize binary orbits over 
time, it is natural to assume that at late times, the inspiralling 
binary is in a "quasicircular" orbit, that is an orbit that is cir- 
cular, apart from the slow shrinkage due to gravitational-wave 
damping. Using numerical relativity it was possible to solve 
the so-called initial value equations of Einstein's theory for 
such quasicircular orbits for a variety of systems, including 
double black holes, double neutron stars, or mixed systems, 
with and without spin. These solutions yielded a set of accu- 
rate values for the orbital energy E and angular momentum 
J as a function of the orbital angular velocity Q.. In Newto- 
nian gravity (which would apply to widely separated binaries), 
these variables would be given by £ = -r)m(GmQ.f^^ /2 and 
J — Grjrn^{GrnQ,)~^^^ . But the systems in these simulations 
were highly relativistic, corresponding to Gm/rc? ~ 0.1 and 
w/c ~ 0.3. Using post-Newtonian expressions for E{Q,) and 
J(r2) valid to 3PN order, it was found that the agreement 



between the PN and numerical results was remarkably good, 
at the level of several percent in most cases [44i |45j. In fact 
it was suggested that some of the systematic differences be- 
tween the PN and the numerical results could be explained if 
the numerical initial configurations actually corresponded to 
slightly eccentric orbits [451 1461 Wf\ . This would change the 
relation between E, J and Q.. Since one is in a highly relativis- 
tic regime, it is not obvious in solving Einstein's initial value 
equations numerically how to choose the initial separation and 
initial angular velocity so as to guarantee a circular orbit ini- 
tially. Subsequently it was discovered that the numerical evo- 
lution of such initial orbits forward in time did indeed display 
small amounts of eccentricity. Numerical methods were then 
developed to fine-tune the initial configurations to ensure the 
desired amount of initial eccentricity. 

In a related development, Favata recently pointed out the 
remarkable effectiveness of the PN approximation for deter- 
mining the final stable circular orbit of two black holes [48] . 

Concluding remarks 

Wigner remarked that the effectiveness of mathematics in the 
natural sciences was mysterious. The unreasonable effective- 
ness of the post-Newtonian approximation in gravitational 
physics is no less mysterious. There is no obvious reason to ex- 
pect PN theory to account so well for the late stage of inspiral 
and merger of two black holes. The Strong Equivalence Prin- 
ciple of general relativity undoubtedly plays a role, by making 
the internal structure of the bodies irrelevant until they begin 
to distort one another tidally. But it does not explain why 
PN waveforms should agree so well with numerical waveforms 
when the orbital velocities are almost half the speed of light, 
or why recoil velocities calculated using PN methods should 
agree so well with those from numerical methods. Our col- 
league Robert V. Wagoner once speculated during the 1970s 
that, because the gravitational redshift effect makes processes 
near black holes appear slower and "weaker" from the point 
of view of external observers, the PN approximation should 
somehow work better than expected, even under such extreme 
conditions. Because of the redshift effect, "strong" gravity is 
not as strong as one might think. But nobody has been able 
to translate Wagoner's musing into anything quantitative or 
predictive. And yet the unreasonable effectiveness of post- 
Newtonian theory will likely be an important factor in the 
anticipated first detection of gravitational waves. 
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Table 1. Current Limits on the PPN Parameters 



Parameter 


Limit 
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2.3 X 10 
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Cassini spacecraft tracking 




4 X 10" 
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VLBI radio deflection 
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3 X 10" 
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perihelion of Mercury 




2.3 X 10 


-4 


no Nordtvedt effect 




10^8 




no anomalous Earth tides 




10-'' 




no anomalies in Lunar, binary 
pulsar orbits 


Q2 


4 X 10" 
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alignment of Sun and ecliptic 


as 


2 X IQ- 


20 


no pulsar "self" accelerations 


Ci 


2 X 10" 


2 


combined PPN bounds 


C2 


4 X 10" 
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no binary "self" -accelerations 


C3 


10"* 




no Lunar "self" -acceleration 


C4 






not independent 
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Fig. 1. Concordance between observations and the post- Newtonian predictions of general 
relativity for the binary pulsar PSR 1913+16. Inset shows the full mi — 7712 plane and the 
intersection region 'a'. The width of each band reflects the error in measuring each parameter. 




Fig. 2. Kick velocity vs. reduced mass parameter Tj for non-spinning black hole mergers. The 
kick vanishes for equal masses (77 = 0.25). Dashed lines indicate estimated uncertainties in 
numerical relativity kicks, while error bars indicate estimated uncertainties in the PN-ringdown 
kicks. 
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